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Abstract 

Chemotaxis is the process by which cells behave in a way that follows 
Cn| the chemical gradient. Applications to bacteria growth, tissue inflammation, 

and vascular tumors provide a focus on optimization strategies. Experiments 
can characterize the form of possible chemotactic sensitivities. This paper 
C ■ addresses the recovery of the chemotactic sensitivity from these experiments 

while allowing for nonlinear dependence of the parameter on the state vari- 
J> ables. The existence of solutions to the forward problem is analyzed. The 

*^ identification of a chemotactic parameter is determined by inverse problem 

techniques. Tikhonov regularization is investigated and appropriate conver- 
gence results are obtained. Numerical results of concentration dependent 
chemotactic terms are explored. 
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1 Introduction 

Biological and ecological research has investigated cell migration. To model cell 
migration, studies have been composed to include migration, diffusion, haptotaxis. 
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and chemotaxis, ||30l|3Tll32lll3llllllll3l. In this paper, the focus is chemotaxis. 
Chemotaxis describes the movement of an organism and/or groups of cells that 
either move toward or away from a chemical or sensory stimulus. In the early 
work by Keller and Segel [30|, chemotactic responses of amoebae to bacteria is 
studied in a cellular slime mold. Bacterial chemotaxis, which describes the ability 
of bacteria to move toward increased or decreased concentrations of attractants 
is analyzed at the macroscopic level through a microscopic model of individual 
cells, |[T8ll47l . It was first observed by Engelmann [14| in 1881. For example, if 
Salmonella typhimurium, a strain of salmonella associated with meat and poultry 
products, is introduced to a petri dish filled with a nutrient, the bacteria will migrate 
outward, consuming the nutrient. As they consume the nutrient, they secrete a 
chemoattractant. After several days, the bacteria will have clustered in the areas of 
high chemical concentration. A structure of concentric rings is usually observed 
experimentally. Chet and Mitchell's work ||8l describes patterns formed from E. 
coli movement toward amino acids. Allweis et. al fT] investigate Vibrio cholerae 
which are inhibited by a pepsin digest that reduces the possibility of the vibrios 
attaching to the intestinal wall. Other authors |[T9l |3l have addressed chemotaxis 
in immune cell motility which when combined with tumor morphology is hoped 
to provide new avenues of treatment strategies. In addition, authors have analyzed 
chemotactic responses in ecology [35] and have investigated mathematical issues 
for the existence of global solutions in multiple dimensions, [34 , 24 , 26] [lOl |9l [B 

Chemotaxis also arises in a variety of medical applications. In particular, it has 
been studied in connection with myxobacteria ||49ll50l . leukocyte mobility in tissue 
inflammation the migration of tumor cells towards bone P4il . and other issues 
in morphogenesis |[39l . Another interesting problem involves the study of vascu- 
lar tumors through angiogenesis. Angiogenesis involves the formation of capillary 
networks of blood vessels that are vital for the growth of tumors. Mathematical 
modehng of angiogenesis ||llSl|7l[IIl|27l|38llSl|42lill3Elhas given new 
insight into tumor structure. Normal tissue, lymphocytes, and other types of cells 
grow at the tumor site or are recruited through chemotaxis. The need to identify 
the nature of this recruitment is at the heart of this paper. The identification of a 
chemotactic term falls under the umbrella of an inverse problem. In principle, we 
can measure certain characteristics of the tumor concentration and use mathemat- 
ical techniques to recover the chemotactic term, in particular the chemotactic sen- 
sitivity, that is driving the tumor growth. To our knowledge, this inverse problem 
approach has only been used in the analysis of chemotaxis models by Dolak-StruB 
and Kiigler 111211 under the assumption that the chemical concentration is explicitly 
known. 

Since there are many applications in which chemotaxis arises, there are also 
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different models of the chemotactic effect. There have been many different expres- 
sions proposed that model chemotactic velocity, see Keller and Segel |[32]| . Lapidus 
and Schiller [36], Ford and Lauffenburger [22], and Tyson et al. (51). This velocity 
is used in a bacterial conservation equation in the formulation of a system of par- 
tial differential equations that governs the particular application. The chemotactic 
sensitivity determines the velocity. Our goal is to develop a technique whereby the 
appropriate chemotactic sensitivity model, and hence chemotactic velocity, can be 
determined from available data. 

In particular, we consider a system of partial differential equations that was de- 
veloped by Oster and Murray ||45]| to model the pattern formation of cartilage con- 
densation in a vertebrate limb bud. A similar system was studied by Myerscough et 
al., [.41] . The numerical solution of similar systems were recently studied by Tyson 
et al. 1571 and by Nakaguchi and Yagi [|42l. Work by Fister and McCarthy [20] has 
shown that the system of partial differential equations can in fact be controlled the- 
oretically through the introduction of a mechanism controlling the number of cells 
being generated. Simulations provide optimal drug treatment programs for patients 
to facilitate the rebuilding of cartilage or the reduction of cancerous tumors. The 
chemotactic sensitivity in |[20l was known and the control parameter was a harvest- 
ing term. Our goal in this work is to identify the chemotactic sensitivity. 

The paper is organized into six sections. In section two, the existence of the 
forward problem is proven. In section three, identifiability of the chemotactic sen- 
sitivity is established using the weak formulation of the state problem. In section 
four, Tikhonov regularization is used to approximate the solution through the use 
of minimization arguments. The rate of convergence of the approximate minimizer 
of the chemotactic sensitivity to the true parameter follows next. In section five, 
numerical experiments provide graphical depictions of the accuracy of the recovery 
of the parameter. In section six, conclusion remarks are made. 

2 Forward Problem 

In this model, u{x, t) and c{x, t) represent the concentration of the cells and the 
chemoattractant, respectively. The cells and the chemoattractant are governed by a 
convection-diffusion equation and a reaction-diffusion equation as 

ut = MAu-V ■{x{u,c)uVc) mnx{0,T) (1) 
ct = DAc+^^-^,c 

u{x,0) = uq{x), c(x,0) = co(x) forxGQ 

l; = i = ondnx{0,T) 
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where v is the outward unit normal. M and D represent the diffusion coefficients of 
the cells and the chemoattractant. The Michaelis-Menten term, ^^p^, represents a 
response of the chemoattractant to a maximum carrying capacity or saturation rate, 
assuming b,h > 0. We incorporate a decay term where fi denotes the degradation 
rate. We assume that there is no flux of the concentrations across the boundary, 
and that the initial concentrations for the cells and chemoattractant are uq{x) and 
co{x), respectively. 

Here, xi^J-, c) is the chemotactic sensitivity which monitors the chemical gra- 
dient attraction of the cells. It is this term that we seek to identify. In B31l4Tll20l 
the term x{u, c) is simply a constant. More generally, xi'u, c) is a linear function 
of u in |El|56l|25l|46l, while in (BOl |32l [36l llll EH it is a nonlinear function of c. 
We assume henceforth that the chemotactic sensitivity has the form x(^i c) = a(c) 
and is a bounded function. We restrict our analysis to the dimensionless system 

Ut = MAu - V • (a(c)nVc) in i7 x (0, T) (2) 
ct = DAc+^,-c 

u{x,0) = uq{x), c{x,0) = co{x) for X G 17 

1^ = 1^ = onaOx(0,r). 

We will establish a technique for the identification of a(c) G A where 



A = {ae H\I) : 



da da 



< A'llci -C2||^2(j) 

L2(/) 



Observe that, with the available data, we can only expect to recover a(c) on the 
interval / = [cmin, Cmax] • The Lipschitz condition on the derivative of the chemo- 
tactic parameter is quite reasonable, since the chemotactic parameter has a rate of 
change that is bounded for bacteria growth, [22] . 

In order to prove identifiability and to establish the rate of convergence ot our 
method, we will need to establish existence of a solution of Q. Using the standard 
notation if '^(17) to represent the Sobelev space W^''^{n), let H^+\9) denote the 
intermediate space between H^{^) and H^+^{9.) for any < < 1. Let D be 
an interval in [0, oo). The space LP{D; X) is the space of measurable functions 
in D with values in the Banach space X. The space C"^{D; X),m = 0, 1, 2, ... is 
the space of m— times continuously differentiable functions in D with values in X, 
while the space C^{D; X),0 < 9 < lis the space of Holder-continuous functions 
in D with values in X. 

Theorem 2.1 Ifuo,co G H^+%n)forO < e < 1, anduo{x) > 0, co(x) > cq > 
on r2, then a real unique local solution u,c of n2\ exists on an interval [0, T] such 



4 



that 

u,ceC'^ {[0,oo);H^+''{n)) DC {[0,T);H^{n)) nC^ {[0,T); L'^{n)) 

with < ei < min(e, ^) and < r] < min (^^^, ^—^^)- The solution satisfies 
the lower bounds 

u{x, t) > 0, c{x, t) > coe"* on [0,T]. 

Proof: Let X = L^{n) x L^{n) and Z = H^+%n) x H^+%n). The system 
([2]) can be formulated as an abstract quasilinear equation 

dv 

— + A{v)v = f{v) < t < oo 
■u(O) = Vq 



on the Banach space X. 
Let 

Clearly G Z. 

We define to be the linear operator in X such that 

^ / -M a{c)u \ f Vu\ f M \ f u 



with domain 



r2/o^ rj2r^^ du dc 



D (AM) = iiie HHQ.) X HHVt); — = — = {)ondn\. 
y on on 

Let the vector f{v) be the function 

Mu ' 



f{v) 

u+l 

Since f{v) is Lipschitz, application of Yagi's work ||58l Thm 2.1 and 3.4] yields 
our result. (See Appendix A for statements of Yagi's results.) □ 
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3 Inverse Problem Statement and Identifiability 



In this section, we begin by establishing the identifiabiUty of the parameter a(c) 
from the available data n(a;, andc(x,t) almost every where in = L'^{{0,T), H^{Q)). 
Note that, in order for chemotaxis to be observed biologically, cells must be present 
and a chemical gradient must exist. This means that u{x, t) and Vc(x, t) must be 
nonzero for a measurable subset of O x (0, T). 

We denote by (ua, Ca) and {uh, Cb) the solution pairs of ^ with chemotactic 
sensitivities a(c) and 6(c), respectively. 

Theorem 3.1 Let {ua, Ca) and {uh, Cf,) both be solutions in W x W of the direct 
problem ([2]) corresponding to a{ca) and b{cb). If Ua = Ub and Ca = Cb almost 
everywhere in i7 x [0, T], then a{c) = 6(c). 

Proof: We consider the weak form of the first equation of the direct problem Q 
for {ua, Ca) and {ub, Cb) and subtract them. 

I [ ^{ua-Ub)(l)dtdx + M [ I {V{ua-Ub))V<t)dxdt 
Jo Jn ot Jq Jq 

= - / [a{ca)uaVca-b{cb)ubVcb]V4>dxdt 
JO Jn 



Since Ua{x, t) = Ub{x, t) and Ca{x, t) = Cb{x, t) a.e., this reduces to 

/ / [a(Ca) - 6(Ca)]UaVCaV(/) dx dt = 0. 

Jo Jq 

By definition, (/> is in = L^((0, T); //^(J7)), then we can choose (f>{x,t) = 
Ca{x, t). Hence, 

T 

/ [a{Ca) - b{Ca)]Ua (VCa)^ dx dt = 0. 

Jn 

Our existence result, Theorem |2.1[ says that > 0. We also employ our biological 
assumptions that / and Vca 7^ on a measurable subset of O x (0, T). Thus 

a(c) = 6(c) 

almost everywhere. □ 



6 



4 Output Least Squares and Tikhonov Regularization 

We wish to identify a function a(c) € A from noisy measurements {zu, Zc) of 
{ua, Ca). Recall that 



A 



a £ H^{I) : 



da da 



< K llci — C2I 



L2(/) 



L2(/) 



where! = [Cmin,Cmax] 

We define 



with 



F{a) = {Ua{x,t),Ca{x,t)) 

F:A^WxW. 



(3) 



In the presence of perfect data (z^, z^), we would solve the non-linear ill-posed 
problem 



F(O^) = {Zu,Zc) 



(4) 



where (n^o , c^o) is the solution of the direct problem with a = a'^. To do this using 
Tikhonov regularization would involve approximating the solution by minimizing 



min||F(a) - (z„, 2;^) H^^^^ + a ||a - a* ||^2(n 



(5) 



where q > is a small parameter and a* is an a priori guess of the true solution 
a^ . In real applications, measurement errors mean that exact data is not available. 
Noisy data is assumed to have an error level 5, which means that 

(■T 



U - Z,, 



dt < 5\ 











f 


c-zi 


Jo 





L2(n) 



dt < 5^ 



(6) 



We assume attainability of a true solution, i.e. if Zc) £ W xW there exists 



a" e A such that 



F(a") = (z„,z,). 



(V) 



In the presence of noisy data {z^,z^), the minimizer a„ G ,A of (|5|) minimizes 

2 



Ja{a) 



F{a)-{zlzi) 



T 



U - Z„ 



II * ii2 

T 

<5 



-a o — a 



dt + 

L'^m JO 

|2 



L2(n) 



dt 



(8) 



for appropriate choices of a G ^ and a. 

We begin by establishing the weak-closedness of the map F{a) 
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Theorem 4.1 If an ^ a* ^ A then Ua„ Ua, and Ca„ Ca, in W. 

Proof: Here, we give the outline of the proof and refer the reader to fTP\ for de- 
tails. Using that the solution to the state system (|2]) is unique, one can define 
Ua„ = u{an) and Ca„ = c{an)- A transformation involving e~^* times each com- 
ponent of the solution pair is made with A to be chosen in order to obtain the 
boundedness of the solution in W. The weak definition of the solution associated 
with the transformed Ua„ and Ca„ in equation ([3]) is analyzed via Cauchy's inequal- 
ity and the boundedness of the coefficients. Using the boundedness (independent 
of n) of the solution pairs, subsequences are extracted that converge weakly to 
and c* . Lastly, comparison results are used so that one can pass to the limit in the 
weak formulation of the solution to show that = and c* = Ca* . □ 
Existence of a minimizer now follows from the lower semi-continuity of 
Ja{a). 



Corollary 4.1 For any data (z^, z^) G W x W, a minimizer o/(|Sj) exists. 

Continuous dependence on the data (zf , ) for fixed a, and the convergence of 
toward the true parameter aP as the noise level 5 and the regularization parameter 
a go to zero also follow from standard results ||48]| . 



Corollary 4.2 For fixed a, the minimizers depend continuously on the data (2^^, 2;^) • 
If a {6) satisfies 

a{6)^0, 6^/a{6)^0 as 6 ^ 



then 



lim 



L2(/) 



4.1 Convergence Rate Analysis 



Although we have noted (without proof) the convergence of the minimizer to 
the true parameter a^, the rate of convergence may be arbitrarily slow. We wish 
to determine a source condition that will guarantee a certain rate of convergence. 
Even when our regularization parameter a is compar able to our noise level 6, such 
a source condition will require assumptions involving u and — a*. 

Recall that we seek to solve the nonlinear problem Q, F{a) = {z^, Zc), where 
F{a) = {ua,Ca)- The true solution is a", and a* is an a priori guess. In order 
to apply the theory of Engl, Hanke, Kunisch and Neubauer ifTTl [T5l , we must 
establish the following: 

• F is Frechet differentiable. 
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• F' is Lipschitz with \\F'{a) - F'{b)\\ < 7 ||a - 6|| , 

• there exists w satisfying the source condition a"^ — a* = F^a"^) 



• 7 < 1. 

In practice, although computing F' and {F')* is not difficult, it can be quite tricky 
to establish the Lipschitz condition on F' with our system of coupled nonlinear 
partial differential equations. Instead, our approach involves developing a source 
condition without imposing differentiability constraints on F. Thus we establish 
0{Vd) convergence. This technique is also found in the work of Engl and Kiigler, 

m. 

Theorem 4.2 Suppose that there exists a function w G L^((0, T); //^(il)) satisfy- 
ing 

w(x, 0) = w[x, T) =0, Aw e l2((0, T);L^(Sl)) 
such that for any ^ ^ A 

r-T 



L^I) 



I I ^{Cao)Uc^oVCaO ■ Vw dx dt. 

Jo Jn 



If a ~ 6 then 



L2(Q) 



+ 



dt = 0{S^ 



and 



0{^5). 



Proof: For clarity, we briefly describe the techniques used in this proof. Using 
that a minimizer to Ja{o) exists, we obtain an upper bound in terms of the error 
level 5 and the norm of the difference in the minimizer and optimal a values. We 
then use our source condition with the weak formulation of the cell and chem- 
ical differential equations to obtain a representation of the inner product of the 
appropriate differences of the approximating minimizers. This allows us to bound 
Ja{a)- Specifically, we use Triangle and Young's inequalities to bound the time 
and spatial derivatives of the differences in the state variables. Integration by parts 
and Holder's inequahty enable us to successfully bound the spatial derivatives of 
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the states in terms of the states themselves. Using the assumptions from A and 
choosing e sufficiently small, we can obtain the error of order \/^ with a ~ (5. 

Since is a minimizer of Ja{o.), we have Jai^a) ^ Jaio.'^)- Using our defi- 
nition of noise level ([6]l, we find that 



< 



T 



U„5 - Z,, 



2 
2 

L2(0) 



+ 



+ 



c„o — Zr 



2 



dt + a 



dt + a \ \a^ — a u j2n\ 
L2{n) II 11^ 

|2 



2 

L2(7) 
2 



< 2(5^ + Q||a° - a* 



(9) 



Adding a ||a'' — a^||^2(/) to both sides of the inequality and using inner product 
properties yields 

rT 







+ 



C^6 Z^. 



dt + a 



<2d^ + 2a{a° -a*,a° -ai 



2 



(10) 



Observe that our source condition 



Jn 



^'(cQo)uaoVc„o • Vw dx dt 



with ^' = a" — , together with the weak forms of the cell equation in the forward 
problem ^ for a'^ and is 



a° -a*,a° -ai 



LHi) Jo Jn 

T 



UnO — u„s I w dx dt 



+ 



Jo. 







V u„o — u„s ■ Vw dx dt 



+ 



Jn 



dx dt 
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and ([TO]l becomes 

rT 




+ 



c„s 



L2(C) 
T 



dt + a 



2 



< 26^ + 2a 







n^o — u^s I w dx dt 



-2aM 



T 







+2a 



Jf7 



■Vwdxdt. (11) 



We bound each integral in ( [TT] ) separately using Triangle and Young's inequalities. 
For the estimates of /i, /2 and I3, we refer the reader to the Appendix for some of 
the details of the bounds used. We find that 



a 



< e 



Jo. 

T 



li^o — u^s I w dx dt 



2 (X 

dt+ — 
L2(n) 2e 



2 rT 



^ dt + eS^ 



and 



aM 



fl 

Jo Jn 



V ( — u^d I • Vw dx dt 



a 
2e 



2 rT 



|A«;||i2(j,) dt + eM^ 











f 




Jo 





dt 



where e is an arbitrary parameter resulting from the use of Young's inequality. We 
utilize the assumptions 



da 



da da 



< K \\ci — C2I 



L2(C) 



L2(Q) 



Green's Theorem, the boundary conditions, and Holder's inequality to obtain the 
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estimate, 

l^3| = 

< £ 



a 



T 



JQ, 

■T 



,2 rT 



1 1 1 1 2 '2 I 

I 9 22 iiii2 2' 



• Vw dx dt 

T 



L2(Q) 



dt 



+^ ll^'"llL°°(f7) 
„2 rT 



K- 



Jo 



2 

L2(Q) 



f 

Jo 



u„s z. 



a f 

+ Y J ll^'^lli^ dt + £\\Vu\\loo^n^K^6'^ + £filii^5'^. 



Grouping terms and relabeling constants, we have 



Jo 



L2(Q) 



+ 



< 26'^ + 2Ci£5^ + 2£C2 



f 

Jo 



L2(f7) 
2 



dt + a 



2 

L2(7) 



dt 



+ 2£C3 



a 

£ 







Jo 






T 


ii2 


MI"^tllL2(Q) 



L2(n) 



_^ (lktlli2(Q)+3||Au;||i2(n) + 2||Vu;||i2(t2)) dt. 



If £ is chosen to be sufficiently small, then for the choices a ^ Sv/e obtain 



Jo 



+ 



and 



L2(n) 



2 

L2(/) 



L2(Q) 



dt = 0{d^ 



2 

L2(0) 



□ 



5 Numerical results 

In order to demonstrate the effectiveness of Tikhonov regularization for this appU- 
cation, we consider several examples. 
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All computations were carried out in MATLAB. The Tikhonov functional 

Ja{a) = j j {\u- zl\^ + \c- zlf^ dxdt + a\\a-a*\\\2(^i) (12) 

was minimized using Isqnonlin, aMATLAB implementation of the Levenberg- 
Marquardt method with line search ||37ll40l . Although it was not tractable to do so 
in the convergence analysis, a gradient based algorithm is appropriate here because 
computing the gradient and its adjoint is straightforward. 

We restrict our discussion to 17 = [0, 1]. Recall that and represent noisy 
data and a* represents an a priori guess of the chemotactic sensitivity a. Cell and 
chemoattractant concentration data on 17 = [0, 1] was generated using pdepe 
with high accuracy. During the computation of Ja(a), cell and chemoattractant 
concentrations u{x,t) and c{x,t) associated with a particular a were computed 
using pdepe with moderate accuracy over coarser space and time meshes than 
those used to simulate data. 

Since Isqnonlin requires objective functions of the form 

l\\F\\l = lY.f"(-)' 
k 

we approximated the first two terms of Ja(a) by 
M r N N •\ 

^ < ^ [u{xi,tj) - Zu{Xi, tj)f + ^ [c{Xi, tj) - Zc{Xi,tj)f > (Ax) (At) 
j=l U=l i=l J 



where Xi = i{Ax) for z = 0, . . . , with Ax = and tj = j{At) for j = 
0, ...,Mwith At = cr/M. 
We approximate a(c) by 

L 

a(c) = '^ak<pk{c) (13) 
fc=i 

where (pk are the usual piecewise linear hat functions defined over a partition of 
[cmin, Cmax]- Notc that any function a(c) can be represented by its corresponding 
vector a. Since the values of Cmin and Cmax may vary considerably for each a used 
during the optimization, we choose instead an interval that is sufficiently large to 
include the range of c for each a considered by the algorithm. In practice, this 
means making a guess, and expanding the interval when c leaves our chosen inter- 
val. 



13 



The penalty term a ||a — a* ||_4 can be replaced by 

L 



^(ai - aX)(\)i 



1=1 



(a-a*)^S(a-a*) 



A 



where the components of the matrix B are given by Bij = {4'i,4'j)A- 

Various strategies for the choice of regularization parameters are discussed in 

II57I. In each of the following examples, we used an L-curve method to choose an 

optimal regularization parameter a. 
Recall that our chemotaxis system is 



ut = MAu - V • {a{c)uVc) in^x (0, T) 

u{x,0) = uo{x), c{x,0) = co{x) for X G O 

ondnx{0,T). 



(14) 



du 
dv 



dc 
du 



A similar system was used by Myerscough et al. f4\] in their numerical simulations 
of chemotaxis in limb-bud development with parameters 

M = 0.25, D = l, a{c) = 2, h = l,b = fi,uo = 1 + e{x),co = 0.5, n = [0, 1] 

where /i G [0, 3000] and e{x) was a bounded perturbation function. In Examples 
1-3, we used the Myerscough parameters with 



e(x) = e-55{^-o.5)2 ^ ^ 0.25, b = ^ = 50. 



Example 1 Consider the chemotactic coefficient used by Myerscough et al. ||4T]| 

a(c, u, X, t) = 2. 

The cell and chemoattractant concentrations associated with this a are shown over 
the time interval [0,0.25] in Figure [T] An initial guess of a = 1 was used. The 
a priori guess was also chosen to be a* = 1 . The parameter a was recovered to 
within 1.461 x 10~^ of the true value at T = 0.25. Commented these next two 
examples out! 

Example 2: Keller-Segel model We consider the nonlinear chemotactic coeffi- 
cient 

a(c) = 2/c 
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proposed in the original Keller-Segel model for chemotaxis ll30ll . The cell and 
chemoattractant concentrations associated with this a are shown at steady state in 
Figure [2] 

From the data, we find that [cmin, Cmax 

] = [0.1794,0.6398] when a (c) = 2/c. 
Since the optimization algorithm will use approximations of other chemotactic 
functions, we attempt to reconstruct a(c) over a larger interval. We found the inter- 
val [0.1, 0.7] to be sufficiently large to include the range of c for each a considered 
by the algorithm. 

An initial guess of a = 15(1 — c)^ was used. The a priori guess was also 
chosen to be a* = 15(1 — c)^. Figure [5] shows the chemotactic function a and its 
recovery of anoise with and without regularization. A regularization parameter of 
a = 10~^ was chosen by an L-curve method. Notice that the regularized recovery 
is quite reasonable over the interval [cmin,Cmax] = [0.1794,0.6398] and that its 
quality degrades, as expected, outside this interval. 



Comments The number of iterations used by the algorithm is quite sensitive to 
the choice of initial function oq and the number of piecewise linear basis func- 
tions used in equation ( 13 1. For experimental data, we must acknowledge that the 
quality of the recovery degrades with increased noise in the data. In certain appli- 
cations such as pattern formation in Escherichia coli or Salmonella typhimurium, 
see Tyson et al. ||5TI . the size of the interval [cmim Cmax] is sometimes too small to 
give adequate information for the recovery of the chemotactic coefficient. This can 
be avoided by taking a larger time interval [0,T]. In numerical simulations, this 
requires a careful choice of numerical method for the solution of the chemotaxis 
system, see [52,1 . An alternative approach is to restrict our measurements to a par- 
ticular time, rather than an interval of time. The efficacy of this approach will be 
discussed in a future paper. 



6 Conclusions 

In this work, we have explored a particular mathematical aspect of the chemotac- 
tic sensitivity within the gradient. The identification of a chemotactic sensitivity 
with functional dependence has been determined. The interesting aspect of this 
work is that, to our understanding, no one has been able to capture the chemotac- 
tic sensitivity information from limited data with dependence on the chemical in 
a system. We have proven the existence of the state solutions in specific Sobolev 
spaces and formulated an inverse problem. We have employed Tikhonov regular- 
ization to recover the chemotactic sensitivity from noisy measurements. In doing 
so, a minimization problem is formed and the necessary convergence results for an 
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approximating minimizer to the true parameter are discussed. 

Another significant result is that we have established a source condition that 
guarantees a particular rate of convergence by imposing a Lipschitz condition on 
the derivative of the chemotactic sensitivity. In practice, this is biologically rea- 
sonable, since the chemotactic sensitivity has a rate of change that is bounded for 
bacterial growth, 1221 . 

Numerically, we have utilized models from Myerscough et. al, f4V\ and Keller 
and Segel, \30l for the studies of the comparison of our proposed work to the 
actual scenarios. With the use of Tikhonov regularization, we have been able to 
recover the chemotactic sensitivity with reasonable accuracy. A biological benefit 
of this knowledge is the ability for one to understand the growth associated with 
chemotaxis within tumor studies, leukocyte dynamics, and bacterial patterns based 
on the specific gradient information that can be recovered from imperfect data. 
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A Appendix 

A.l Yagi's existence theory 



For more thorough understanding of Theorem 2. 1 we include the parabolic prob- 
lem convention used in Yagi, ll58l . 

ut = V ■ {a{u,p)Vu) — uh{p)Vp) ini7x(0,oo), (15) 

Pt = dAp + uf{p) — g{jp)p in X (0, oo), 

0) = no(x), p{x,{)) = pq{x) forxerj 

i; = §2 = onai7x(0,oo). 

In Yagi's theorem ll58l Thm 2.1], it states 
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Theorem A.l Letuo,po £ i7 ^+^0 /or < eo < 1> and uo{x) > 0, po{x) >> 



^0 > o?i 0,. Assume a real local solution {u, p) to {15) exists on the interval [0, S\ 
such that 

u,pec ([0, s);H^-^''{n)) n c ([0, s);H^{n)) n ([o, sy, L\n)) 

for some E\ > 0. In addition, assume that p satisfies p{x, t) > onO, x [0, 5] and 
an estimate 

\\v{t)\\m < At^^onO <t<S, 
for some 82 > and constant A. Then 

u{x,t)>0, p{x,t)>p{t) forall{x,t)£Tix[0,S], 

where p denotes a positive function defined as the global solution to an ordinary 
differential equation, 

dp 

= —g{p)p onO < t < 00 (16) 
p{0) = 60 >0. 



For further connection to our work, we utlilize the continuation to a unique 
solution that Yagi developed |58, Thm 3.4] with 

Theorem A.2 Let uq.pq £ H^+^"{n)forO < eo < 1, and uq{x) > 0, ^0(2;) > 
^0 > on $7 and letO < r] < 13— a. Then, in the function space, C ([0, 00) ; H^'^^^ (J7 



the problem ( 15 ) possess a unique local solution 

u,pec ([0, s);H^{n)) n ([0, sy, L^{n)) 

with the lower bounds 

u{t)>0, p{t)>p{t) forte[0,S], 



where p{-) is a positive function defined by equations in ( 16 1. The interval [0, S] on 
which the solution exists at least is determined by the norms 1 1 tio 1 1 ji/i+eo > I |Po 1 1 ^ i+^o 
and by the initial lower bound 60. 

It is to be noted that by this theorem from Yagi's work that a maximal solution to 



( 15 1 can be uniquely defined in the space C ([0, 5); if^+^i (17)) for < r/ < /3-a 
for each uo,po such that uq,pq G //^"'"^^(O) for < eo < 1> and uo{x) > 0, 
Po{x) > Co > (5o > on rj. 
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U{X,t) C(X,t) 




Figure 1: Cell and Chemical Concentrations over time interval [0,0.25] with 

o(c, u, X, t) = 2. 



u(<,t) c(x,t) 




Figure 2: Cell and Chemical Concentrations over time interval [0, 0.25] with 

a(c, u,x,t) = 2/c 
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a(c)=2/c, 20 points, alpha=1e-5 




Figure 3: True chemotactic coefficient a{x) = 2/c and its recovery with no regu- 
larization (left) and with a regularization parameter of a = 10~^. 
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